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ABSTRACT 


We investigate small-scale dynamo action in the solar convection zone through a 
series of high resolution MHD simulations in a local Cartesian domain with 1 R & (solar 
radius) of horizontal extent and a radial extent from 0.715 to 0.96 R Q . The dependence 
of the solution on resolution and diffusivity is studied. For a grid spacing of less than 
350 km, the root mean square magnetic field strength near the base of the convection 
zone reaches 95% of the equipartition field strength (i.e. magnetic and kinetic energy are 
comparable). For these solutions the Lorentz force feedback on the convection velocity is 
found to be significant. The velocity near the base of the convection zone is reduced to 50% 
of the hydrodynamic one. In spite of a significant decrease of the convection velocity, the 
reduction in the enthalpy flux is relatively small, since the magnetic held also suppresses the 
horizontal mixing of the entropy between up- and downflow regions. This effect increases 
the amplitude of the entropy perturbation and makes convective energy transport more 
efficient. We discuss potential implications of these results for solar global convection and 
dynamo simulations. 
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Introduction 


Turbulent thermal convection fills the solar convection zone due to its superadiabatic 
stratification. Thermal convection under the influence of rotation leads to angular 
momentum transport and maintains large-scale mean flows, such as the differential rotation 
and the meridional flow. These flows in combination with helical turbulence play a 


crucial role for maintaining the Sun’s global magnetic field through a la rge-sca 


e dyn amo, 


Steen 

beck & Krause 

1969a 

b; 

Krause & Radler 

1980; 

Osscndriiver 

2003; 

Charbonneau 


2005; 

Miesch 

2005 

). In addition, also non-helical turbulent thermal convection itsel 

' has the 


ability to maintain a turbulent magnetic field (small-scale dynamo e.g. 


Batchelor 


19501) ) 


Small-scale dynamos maintain magnetic field s on scales comparable to or smalle r than the 


energy carrying scale of turbulence (see, e.g., 


Brandenburg fe Subramanian 


and are sometimes also called “fluctuation dynamos” (see also 


,2005 


, s ectio n 5) 


Brandenburg et al 


Numerical simulations of a small-scale dynamo were presented first by 


2012 


Meneguzzi et al 


(119811 1 for non-stratified forced turbulence, followed by high resolution small-scale dynamo 


simu 

(e-g- 


ations in the context of galaxies, stars. 

he Sun as well 

as idealized forcing 

urbu 

Brandenburg et al. 

1996; 

Hawlev et al. 

1996; 

Cattaneo 

1999; 

Haugen et al. 

2004 


Schekochihin et al. 


ence 


20041 ). Small-scale dynamos are excited when the magnetic Reynolds 


number exceeds a critical value, typic ally near 100, which is found to be dependent on th e 


magnetic Prandtl number Pm = v In (IHaugen et ah 


Iskakov et al. 


2007 : 


Schekochihin et al 


2004 : 


Schekochihin et al. 


2004 , 


2005 


20071 ). The critical magnetic Reynolds number is 


found to be larger by a factor of a few for the P m -C 1 regime that is most relevant for 
stellar and planetary dynamos. Compared to large-scale dynamos, small-scale dynamos 
have in their kinematic phase small growth time scales that are determined by the eddy 
turn-over time scale near the resistive scale (Pm -C 1) or viscous scale (Pm 3> 1), i.e. the 
growth rate is strongly dependend on the magnetic (and viscous) diffusivity. While the 
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onset and kinematic growth of small-scale dynamos has been studied in great detail, the 
non-linear saturation p hase in pa rticula r for the Pm <C 1 regime has been addressed only 


very recently (Brandenburg 


2011 


201411 


The non-linear saturation regime of a small-scale dynamo is most relevant for a 
theoretical explanation of the observed mixed polarity held in the solar photosphere 
outside active regions, the so-called “quiet sun”. Small-scale dynamo simulations that 
include more realistic effects, such as partial ionization and radiation in the uppermost 


few M m of the convection zone were p resented recently by 


Pietarila Graham et al 


12010) 


Vogler fc Schiisslerl (12007(); 


ReinpeJ (2014) found that explaining the observed level of 


quiet Sun magnetic held requires solutions that reach a magnetic held strength close to 
equipartition only a few Mm beneath the photosphere. Small-scale magnetic held with such 
strength has potentially a significant influence on convective dynamics. 

The above mentioned photospheric dynamo simulations cover only the uppermost few 
Mm of the solar convection zone and typically use open bottom boundary conditions to 
“mimic” the influence from the deeper convection zone. Due to the use of open boundary 
conditions these simulations cannot determine the saturation held strength of a small-scale 
dynamo in the photosphere without assumptions about the magnetic held strength and 
structure in the deeper convection zone. We p resent here small scale dynamo simulations 


Rempelj (2 0141 ). but focus on the lower convection 


that are very similar to those presented by 
zone ranging from 0.715 to 0.96 Rq. 

The purpose of this study is twofold: Firstly we want to investigate whether the 
assu mption about the magnetic held strength in the convection zone made in 


Rempel 


(120141) in order to explain the observed quiet Sun held strength are realistic. Secondly we 


want to investigate how a strong (about equipartition) small-scale magnetic held influences 
convection and convective energy transport throughout the convection zone. 
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In global dynamo simulations small-scale and large-scale dynamo action is present at 
the same time, however most global dynamo simulations have insufficient resolution to 
properly address the small-scale dynamo: The sun has a circumference of 4400 Mm at the 
surface, which is large compared to the typical scales of the convective structures and the 
local pressure scale height (~ 60 Mm at the base of the convection zone ). Most investigations 


of the global solar dynamo have a 


Augustson et al. 

2013; 

Fan & Fang 

2014) 


grid spacing larger than 4 Mm flBrun fe Toomre 


2002 


2 0141 ). This resolution is not sufficient to resolve the 


inertial range of the turbulence, which is important for the small-scale dynamo. While the 
small-scale dynamo is not completely absent, its efficiency and saturation field strength are 
likely significantly underestimated. 


Recently 


Hottaetah 


(120141 ) (hereafter Paper I) studied a small-scale dynamo in 


a global setup without rotation (in order to rule out large-scale dynamo action by 
construction) and achieved a grid spacing of 1 Mm and found rather efficient small-scale 
dynamo action throughout the convection zone. Magnetic field of 0.15R e q — 0.25R C q is 
maintained by the small-scale dynamo, where B cq = y/4:7igv ims is the equipartition magnetic 
field strength. Although they could see feedback from the generated magnetic field on the 
velocity, the effect is not significant for that resolution. 

The purpose of this study is to increase the resolution further by limiting the 
size of the calculation domain and to investigate the efficiency of the small-scale 
dynamo and its feedback on convection using a grid spacing as small as 350 km (a 


global convection simulation with this resolution would require a g ric 


600(r) x 6000(0) x 12000(0)). In addition we use (similar to Rempell (2014)) a LES (Large 


of approximately 


Eddy Simulation) approach, which drastically reduces diffusivities on the resolved scales 
of the simulations. Our approach implies that the effective magnetic Prandtl number 
is close to unity and we do not carry out a parameter survey for the magnetic Prandtl 
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number in this study. Our focus is on evaluating the potential influence of an efficient 
small-scale dynamo on convective dynamics, which is currently not resolvable in global 


models. This is q uite differen 


from studies of glo 


jal dynamos in the Sun and in rapidly 


rotating systems ( 

Brun et al. 

2004; 

Ghizaru et al. 

20 

10; 

Brown et al. 

20 

11; 

Fan & Fang 

2014; 

Christensen et al 

1999 

Busse & Simitev 

2011; 

Yadav et al. 

2013; 

Schrinner et al. 

2014; 

Yadav et al. 

2014 

, where a significant fraction of the magnetic energy resides on t 


largest scales. 


Model 


We solve the three-dimensional m agneto 


speed of sound technique (RSST: 


lydroc 


Hotta et al. 


ynamic equations with the reduced 


2012b ) in Cartesian geometry (x,y,z). In our 


definition, ^-direction is the vertical direction. The equations are expressed as 


9pi = —ly-(pov), 


Po 


dt 

dv 

~dt 

OB 

~dt 


B 2 \ 


= Po(v ■ V)v - V {pi + —1 + V • 

= V x (v x B), 


BB 

47T 


Pipe x , 


T ® Sl T ( Y7\ -L ^ { dT 0 

P0 T 0~gj- = -PoTflV ' V)Si + — ( K r-P0Cp — 


+ r, 


'dp\ ,fdp 


Si. 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 


The subscripts 0 and 1 denote that background and time-dependent perturbed values. 
Except for the geometry, the setting is similar to that presented in Paper I. While a 
limited horizontal extent of the calculation domain is required in order to increase the 
resolution, we also intend to avoid boundary effect from the side boundaries. Thus we use 
a periodic boundary and Cartesian geometry. One disadvantage of this setting is that it 
leads to the incorrect energy flux in the upper part of the domain. The solar energy flux 





















































at the base of the convection zone (r = 0.715i? o , where R Q is the solar radius), is about 
1.21 x 10 11 erg s -1 cm -2 . This is about twice as large as the flux found in the photosphere 
due to the expansion factor of r 2 . We should keep in mind that the obtained convective 
velocity will be larger by less than factor of \/2 in the upper part of the domain. The s caling 


law of v 2 s oc F is found in some studie s even in rapid rotation situation (IChristensen 


2010 


Schrinner et al. 


2014 


Hotta et ah 


2015J), which is not far from mixing length theory, which 


predicts the dependence of v 3 oc F (e.g. 


Stix 


2004j). The background stratification, p 0 , p 0 , 


and T 0 , is obtained by solving the hydrostatic equation in the adiabatic atmosphere. The 
gr avitational acceleration q and the radiative diffusivity K r are adopted from the Model 


S ( Chr i stensen-Dalsgaard et al . 


19961 ). The obtained stratification is almost same as the 


Model S (see Fig. 1 of Paper 1). The same form for the surface cooling T is adopted as in 
Paper I. The thickness of the cooling layer is two pressure scale heights at the top boundary, 
which is 18.8 Mm. We use a realistic equat ion of state, which include the partial ionization 


effects using OPAL repository ((Rogers et al 


19961 ) is used in eq. ([5]), the plasma is almost 


completely ionized in the calculation domain of this study. Thus the equation of state is 
almost the same as that for an ideal gas. The factor for the RSST is set for keeping the 
reduced adiabatic speed of sound c s = c s /£ constant as 

c s (x) 


f(z) = £o- 


( 6 ) 


Cs ('C'min) 

where c s (x) = \J(dp/dp) s and x min the original adiabatic speed of sound and the location 
of the bottom boundary. We use £ 0 = 150 and the reduced speed of sound c s is 1.5 km s _1 . 
Considering the validity of the RSST, i.e. u r ms/c s < 0.7, we can properly treat the convection 
of u rms < 1.05 kms -1 , where u rms is root mean square (RMS) velocity. In addition to 
in creas ing the effective resolution, we also adopt a less diffusive artificial viscosity suggested 


m 


Rem pell (120141 ). We note that this modification has a significant effect on resolving the 


small-scale magnetic held and leads to significantly larger kinematic growth rates, i.e. more 
efficient small-scale dynamos. 































Our calculation domain extends (0.715i2 o , 0, 0) < ( x,y,z ) < (0.96 R Q , R Q , R & ). The 
periodic boundary condition is adopted for the y- and ^-directions (x is vertical). The 
symmetric boundary condition is adopted for density and the entropy perturbation at the 
top and bottom boundaries. Impenetrable and stress-free boundary condition are adopted 
for the velocity {v x = dv y /dx = dv z /dx = 0). At the top boundary, the horizontal magnetic 
field is zero (B y — B z — 0). At the bottom boundary a symmetric boundary condition is 
adopted for the horizontal magnetic field. Then the vertical magnetic field B x is calculated 
to satisfy the V ■ B = 0 constraint at top and bottom boundaries. In the calculation 


domain the c 
suggested by 


i vergence free condi tion for the magnetic field is maintained with the method 


Dedner et al. 


( 120021 1 . 


Our investigation is focused on the solar convection zone. Since we investigate the 
contribution from a small-scale dynamo alone, we do not include rotation. In this setup 
the stratification and the luminosi ty are well determined by a solar standard model 


(jC hr is tense n-D alsg aard et al. 


19961 1. The free parameters in this study are domain extent 


and resolution (which translates into diffusivities since we use an LES (Large Eddy 
Simulation) approach). The typical scale of stratified convection is determined by the local 
pressure scale height, which reaches a value of ~ 60 Mm at the base of the convection zone. 
Our horizontal extent of 700 Mm is large enough to expect at best a weak influence. The 
vertical extent of our domain from 0.715 to 0.96 Rq was chosen to capture most of the 
convection zone, but leaves out scales near the top, which are too small to be resolved. 
Since we focus here on small-scale dynamo simulations it is not sufficient to just resolve the 
scale of convection, we need to be able to resolve also a significant fraction of the turbulent 
energy cascade. The effective magnetic and fluid Reynolds numbers are tied to the grid 
resolution in our study due to the use of an L ES approach. Estimating the fluid and 


magnetic Reynolds numbers is a difficult task (iRempe] 


201 if). We present one simulation 


with an explicit magnetic diffusivity for which we find a value of Rm in the range from 300 
















to 1200 (see following paragraph). In our highest resolution LES simulation the small-scale 
dynamo has an about 100 times larger growth rate, which indicates a substantially (at least 
2 orders of magnitude) larger ’’effective” magnetic Reynolds number in those setups. (Note 


that scaling arguments suggest for the growth rate of a small-scal e dynamo y vR ; ! 111 . 


which would imp 


(| 2010 |); 


y an effective Rm a factor of 10 4 larger, however, 


Pietarila Graham et al. 


Rempel 020141 ) found more a 7 ~ Rm relationship for the resolutions currently 


accessible by numerical simulations.) 

We carry out four hydrodynamic cases and corresponding four MHD cases. The 
difference is mainly in the resolution (see Table CO). The characters “H” and “M” denote 
hydrodynamic and MHD. Cases H256D and M256D include explicit thermal diffusivity, 
kinetic viscosity and magnetic diffusivity in order to compare the result with previous 
solar global dynamo calculations. The thermal diffusivity( k), kinetic viscosity(^) and the 
magnetic diffusivity( 77 ) are added as: 

<9v 


P 0 


P0T0 


dt 

<9B 

dt 

dsi 

~dt 


-V D, 

- V x (??V x B), 


V • (kpoTqVsi) 


+2pqu 


&ij &ij 


O (V • v) 


+ ^ VXB) 


(7) 

( 8 ) 

(9) 


where D is the viscous stress tensor defined as 


Dij = -2p 0 u 




j(V ■ v)<% 


( 10 ) 


Fan fe Fan g 


and C jj is the strain rate tensor. The distribution of v and 77 are similar to 
( 20141 ). At the top boundary, the value is 1 x 10 12 cm 2 s _1 and these decrease with depth 


following a 


U/Pn profi le. This 


calculation (Miesch et al. 


2000 


setting is common in ASH (Anelastic Spherical Harmonic) 


Brun et al 


20041 ). In this study we adopt the same value 
for the thermal diffusivity as the other values (k = lx 10 12 cm 2 s _1 ). We use the definition 
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of the magnetic Reynolds number Re = v ims d c /r] (e.g. 


Gastine et al. 


20121 ). where d c is 


the thickness of the convection zone. The estimated magnetic Reynolds number is 336 
and 1250 at the top and bottom of the calculation domain, respectively. u rms is obtained 

le com parabl e valu e to the recen t high reso 


from the resu 

t. This is t 

(Brown et al. 

2010 , 

r-H 

i—1 

O 

CM 


Gastine et al. 


2012 ; 


Fan fc Fang 


2014 


reso 

ution calculations 

Jones 

2014; 

Yadav et al. 


2 0141 ) All calculations started from a hydrodynamic case and were evolved for 100 days. 


Then, weak random magnetic field ( B z ) with amplitude of 100 G is added to initiate a 
small-scale dynamo. The imposed magnetic held is uniform in the direction and random 
in the x — y plane. The name of cases becomes “M” correspondingly. In the highest 
resolution case (H2048 and M2048), the grid spacing is smaller than 350 km. 


Results 


3.1. Structure of velocity and magnetic field 


We start our discussion with cases H256D and M256D. These mimic the currently 
available global dynamo calculations. The grid spacing is 2700 km and the explicit 
diffusivities are included. We note that almost isotoropic grid spacing is used in this 
paper, i.e, Ax ~ Ay = Az. We calculate case M256D for 2000 days in order to cover 
its long time scale to the saturated phase. Figs. [T^ and b show the contour of vertical 
velocity v x and vertical magnetic held B :r at r = O.95i? 0 , respectively. As is often seen in 


global dynamo calculations (e.g. 


Brun et al 


20041 ). broad uphows are surrounded by thin 


downflows. The vertical magnetic held is concentrated in the downflow region. Fig. [2] 
shows the temporal evolution of magnetic energy averaged over the computational domain 
in case M256D. The e-folding time scale for magnetic energy during the kinematic phase 
is 112 days. This is longer than the convective time scale at the base of the convection 
zone, which is H p /v rras ~ 10 day, where H p and u rms are the density scale height and 
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the RMS velocity, respectively. We note that when we double the diffusivities, i.e., 
k = v = r) = 2 x 10 12 cm 2 s _1 , we do not find an excited small-scale dynamo. Fig. |3Jr shows 
the distribution of the equipartition magnetic held strength B cq = \/Anp 0 v rms (black lines) 
and the RMS magnetic held B rms (red line). The dotted black line shows the result of the 
hydrodynamic case (H256D). The RMS magnetic held is about 3000 G and 500 G around 
the bottom and top boundaries, respectively. The small difference of solid and dotted black 
lines shows that the influence of the Lorentz feedback on the convection is very small in this 
case. The dynamo in case M256D maintains 0.08R eq — 0.20R cq magnetic held. The result 
indicates that the small-scale dynamo in currently available global dynamo calculation is 
not efficient and that the generated magnetic helds have only a negligible effect on the 
convective velocity. These findings concern only to the contribution from small-scale held. 
The large-scale held of a dynamo in ra pidly rotating systems can have a significan tly effect 


on the turbulence in some studies (e.g. 


Christensen et ah 


1999 


Yadav et al. 


2 0141 ). In the 


following investigation we decrease the grid spacing to values of less than 350 km, which 
currently cannot be reached on the global dynamo calculation. 


All simulations with higher resolution do not use any explicit diffusivities. The 
investigation on hydrodynamic cases (H512, H1024, and H2048) is done for the time 
interval 80-100 days. Fig. [4] shows the distribution of the vertical velocity at r = O.95i? 0 
(panels a and b) and r = O.8OR 0 (panels c and d). Compared with case M256D (Fig. Q3i), 
convection cells are distorted by small-scale turbulence (Fig. Hb). Following the increase of 
pressure scale height, the typical convection cell size becomes large in the deeper layer. In 
addition, small-scale turbulence exists much even in the upflow region (Fig. |Ud). Figs. EK, 
b, and c show the spectra of kinetic energy at r = O.95R 0 , O.8i? 0 , and O.72R 0 , respectively 
for cases H256D, H512, H1024, and H2048. The dotted line indicates a Kolmogorov 
slope of power-law index —5/3. In the upper layer (r = O.95i? 0 ), the spectra of kinetic 
energy obey the power-low distribution with the index of -5/3 in cases without explicit 
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diffusivities. In the deeper region, especially in r = O.72i? 0 , there is a depression around 
k/2n = 2 x 10" 1 Mm -1 which is largest in case H2048 (red line in Fig. [5fc). 


Although we calculate the cases M512, and M1024 for 900 days, no significant change 
can be seen after 300 days. We calculate case M2048 for 300 days and compare the results 
at the time of 280-300 days. Fig. O shows the temporal evolution of magnetic energy 


averaged over the whole calculation domain. As stated in 


Rempell (12014T) . the e-folding time 


scale for the growth rate of magnetic energy significantly depends on the resolution. The 
time scales are 3.24, 1.97, and 0.99 days at case M512, M1024, and M2048, respectively. In 
addition, the saturated magnetic energy is one order of magnitude higher than that in case 
M256D. The magnetic energy in cases M1024 and M2048 converges around t = 300 days. 
Figs, d and [8] show the distribution of vertical velocity v x and vertical magnetic field B x at 
r = O.95i? 0 and r = O.8i? 0 , respectively. Compared with the hydrodynamic case (H2048: 
Fig. [4j), the appearance of the convective flow is smoother. The small-scale turbulence in 
upflow regions is strongly suppressed (see Figs. |4jl and Eb) . In the upper layer, the location 
of small-scale features in v x have rough coincidence with the location of strong vertical 
magnetic held B x . Small-scale magnetic held with mixed polarity is preferentially found in 
downflow regions. In deeper layers, the coincidence between downhow and strong magnetic 
held is less pronounced. The magnetic held is distributed rather uniformly. Joint PDFs for 
vertical velocity (v x ) and vertical magnetic held (B x ) also show this feature in Fig. |9j While 
the distribution of magnetic held smoothly connects between up- and downhow region in 
the deeper layer (panel b), a steep transition is seen near the upper boundary (panel a). 
We note that the ratio of the RMS values of vertical magnetic held in uphow to that in 
downhow (~ 0.5) is not different between these layers. 


This finding is conhrmed with the spectra of magnetic energy E mag and kinetic energy 
Akin in Fig. [TOJ The black and blue lines show A kin and A mag in case M2048. The red 
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lines show E^,, in case H2048. Panels a, b, and c show the results at r = O.95i? 0 , O.8i? 0 , 
and O.72i? 0 , respectively. Fig. HUH shows the ratio of E mag to E kin . The black, bine, and 
red lines show the result at r = O.95i? 0 , O.8i? 0 , and O.72i? 0 , respectively. In every layer, 
super-equipartition magnetic held is found at small scales. There is depth dependence of 
the scale where the magnetic energy exceeds the kinetic energy. The scales are 8.1, 15, and 
31 Mm for r = O.95i? 0 , O.8i? 0 , and O.72i? 0 , respectively. In Paper I, it is argued that there 
are two factors that influence the dependence of dynamo efficiency on the depth. One is 
the negative Poynting flux in the convection zone. In the upper convection zone, magnetic 
energy is rapidly transport downward and accumulates in the lower convection zone. This 
makes small-scale dynamo action more challenging near the top and less challenging near 
the bottom of the domain. The other factor is the dependence of the intrinsic spatial scale 
of the convection flow on the depth, i.e., the convection in deeper layer has larger spatial 
scale. For a fixed grid spacing convective structures are better resolved near the bottom 
of the domain. The super-equipartition magnetic energy is 6 times larger than kinetic 
energy in the diffusion scale {k/2n > 0.5 Mm -1 ). This likely indicates a numerical magnetic 
Prandtl number larger than unity in our calculation. In upper layer (r = O.95i? 0 : Fig. 
HOhh small-scale kinetic energy is selectively suppressed by magnetic held and the larger 
scale kinetic energy remains. On the other hand, all the scales are suppressed in the deeper 
layer at r = O.72i? 0 (Fig. flOfc). 

Fig. [TTh shows B cq and B ims . The solid and dash-dotted lines show B eq and £> rms in 
case M512 (black), M1024 (blue), and M2048 (red), respectively. The dotted lines show B cq 
in hydrodynamic cases H512 (black), H1024 (blue), and H2048 (red). Although around the 
base of the convection zone, the strength of the magnetic held between cases M1024 and 
M2048 converges, the reduction of B eq is larger in M2048. The higher resolution simulation 
has more kinetic energy on small scales (recall the depression in Fig. [5h). Since the Lorentz 
feedback has a larger effect at smaller scales, the reduction of B cq is larger in case M2048. 


14 


Fig. ITTb shows the ratio of B vms to B cq in cases M512, M1024, and M2048. At the base of 
the convection zone, case M2048 achieves 95% of the equipartition field strength. 


Fig. IT2l shows the horizontally averaged work densities of pressure/buoyancy (Wd) and 
the Lorentz force (Wi) defined as 

/ v ■ (Vpi + pige x )dS 


W d = 


W\ = 


fdS 

fv -(j x B )dS 

T dS ’ 


( 11 ) 

( 12 ) 


where S is the horizontal plane. The work by pressure and buoyancy almost converges 
between M1024 and M2048. The conversion rate from kinetic energy to magnetic energy 
(W\) is 36%, 61%, and 76% of the pressure/buoyancy work (Wd) in the calculation domain 


in cases M512, M1024, and M2048, respectively. The magnetic Prandtl number is c 


unity in this study since we only use numerical diffusivity. Recently 


Brandenburg (2011 


ose to 


2 0141 ) found that for efficient dynamos with Rrn 1 the work by the Lorentz force and 


consequently the ratio of viscous to resistive energy dissipation is dependent on the magnetic 
Prandtl number. In particular in the regime 1 <C Rm -C Re, which is most relevant for solar 
and stellar convection zones most of the energy dissipation happens through resistivity. 
This means that the value Wi/Wd increases with lowering the magnetic Prandtl number 
and should be close to unity (i.e. the small-scale dynamo is maximally efficient). Note that 
in global dynamo simulations with moderate Rm, where an efficient small-scale dynamo is 


not present, the value Wf/Wd (most 


parameters such as Rossby number (tSchrinner 


y d ue to a large-scale dynamo) can depend on other 


20131). 


Fig. fl3l shows the RMS velocity for all cases. In panel a, the solid (dotted), black, 
blue, and red lines show the results in cases M512 (H512), M1024 (H1024), and M2048 
(H2048). Figs. FTTH ), c, and d show the ratio of RMS velocity in the MHD case to that in 
the hydrodynamic case. The black, blue and red lines show the result in cases H(M)512, 
H(M)1024 and H(M)2048. The panels b, c, and d show the RMS velocities including 
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all three components, the x-component, and the horizontal component (uh = \/v x + Vy), 
respectively. As shown in Fig. [TU Fig. [T3]also shows that the reduction of the velocity with 
magnetic field does not converge with M1024 and M2048. In the highest resolution case 
M2048, the RMS velocity is reduced to 85% and 50% at the top and bottom boundaries. 
Figs. fl3b and d show that horizontal velocity is more suppressed by magnetic field, which is 
a key feature for understanding the energy transport in the following discussion. While the 
suppression of v x is converged, the suppression of u h shows strong resolution dependence. 
As shown in the spectra (Fig. flOb and c), the kinetic energy is reduced independently of the 
spatial scale from the middle to the base of the convection zone. As a result, the unsigned 
mass flux is also suppressed in the deeper region (Fig. [14]). This indicates that the amount 
of over turning is reduced. It is thought that this is caused by the large-scale magnetic field 
component, since the reduction does not depend on the resolution significantly. 


Figs. [T5h . and b show PDFs for v x , and v y , respectively. The dotted and solid lines 
show the results in cases H2048 and M2048. The black, blue and red lines show the result 
at r = O.95i? 0 , 0.8 Rq, and O.72i? 0 , respectively. As shown in Paper I, the vertical velocity 
has a asymmetric distribution and the distribution of the horizontal velocity is almost 
Gaussian. Larger velocity values are more strongly reduced at all depths. Fig. [161 shows 
the kurtosis /C and skewness S , which show intermittency and asymmetry of distribution, 
respectively. As explained above, the horizontal velocity has almost Gaussian distribution 
(/C 3 and S ~ 0) and vertical velocity (v x ) has higher intermittency and asymmetry in 

the hydrodynamic case (H2048: dotted lines). When the magnetic field is included (M2048: 
solid line), intermittency is increased especially at the bottom part of the convection zone. 
This indicates that the higher velocity is selectively suppressed by the magnetic field. Figs. 


IT5fc and d s 
to those in 


l ow P DFs for B x , and B y , respectively. The PDFs of magnetic field are similar 


Rempel (120141 ). Due to the boundary condition at the bottom, the vertical 


magnetic field at r = O.8_R 0 is larger than that in r = O.72R 0 . The peak magnetic field 
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strength is around 65 kG near the base of the convection zone. Although we do not generate 
flux concentrations that are large enough to lead to flux emergence, the ach ieved m agnet ic 
field is sufficiently stro ng in order to reproduce the active region features flWeber et al. 


2011 


Hotta et al. 


2012a) . The vertical field around the top boundary (black line in Fig. 


[T5b ) shows an extended tail. Although this is also seen in photospheric simulations, our 
result is likely strongly influenced by the impenetrable boundary condition at the top. 


3.2. Energy transport in high resolution calculations 

Fig. [171 shows the horizontally averaged vertical energy fluxes. The definition of the 
enthalpy flux F e , the kinetic flux F k , the radiative flux F r , the Poynting flux F m , and total 
flux F t are shown as 


F e 

f , PoPi \ 

= Po e i + Pi 

\ Po ) 

v r , 

(13) 

F k 

1 2 

= g p ° v Vr ’ 


(14) 

F r 

dT 0 r 

— fiypoCp ^ F / 

^ *^min 

T (x)dx, 

(15) 

F m 

= 7 '-{EyB z - E Z By), 

47T 


(16) 

F t 

— F e + F k + F r + F m , 


(17) 


where e\ is the internal energy and E = — (v x B)/c is the electric field. The artificial 
surface cooling is included in the radiative flux F v in order to avoid complexity in the 
figure. We find a downward directed Poynting flux in our study that reaches 8% of the 
solar energy flux. This is three times larger than the value found in Paper I. Whereas the 
reduction in the RMS velocity is large (25% reduction in the middle of the convection zone), 
the reduction of the enthalpy flux is relatively small (12% at maximum), which implies 
an increase of the entropy difference between up- and downflows. Fig. [T8l shows PDFs of 
the entropy perturbation. The dotted and solid lines show the results in cases H2048 and 
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M2048, respectively. The black, blue and red lines shows the results at r = O.95i? 0 , O.8i? 0 , 
and O.72i? 0 , respectively. In every layer, the entropy perturbation, especially in negative 
part, has larger amplitude in case M2048. Figs. [TUk and b show joint PDFs with vertical 
velocity v x and the entropy perturbation .si in cases H2048 and M2048, respectively. In the 
MHD case (Fig. fl9b). the larger amplitude of the entropy perturbation is achieved in the 
downflow region. This improves the overall efficiency of the convective energy transport. 
Figs. I20h and b show the contour of the entropy at r = O.8i? 0 . Unlike the vertical velocity 
(Figs. 0k and [8k), the presence of small-scale magnetic field leads to smaller features in the 
entropy perturbation. In addition the area of negative entropy perturbation is reduced by 
20% at maximum. The area of negative entropy perturbation is 30% and 23% at the middle 
of the convection zone in the cases H2048 and M2048, respectively. This is also seen in 
the spectra of entropy in Fig. 0U In the deeper layer in particular the small-scale entropy 
perturbation is amplified. As a result, the enthalpy flux in the middle of the convection 
is shifted to smaller scales (Fig. [22]) . In other words, strong magnetic field enhances the 
small-scale energy transport in the middle of the convection zone. Fig. [23] shows the RMS 
values of the entropy in cases H2048 and M2048. In the downflow region (blue lines), the 
amplification of the entropy perturbation is more pronounced. The effect of the magnetic 
field in upflow region (red lines) is mostly seen around the base of the convection zone. 
Fig. 1241 shows the kurtosis and skewness of the entropy in cases H2048 and M2048. The 
intermittency and asymmetry of the entropy are increased by the magnetic field. Figs. 
[T3b and d show that the horizontal velocity is more suppressed by the magnetic field than 
the vertical velocity. In addition, Fig. [25] shows the ratio of horizontally averaged Isii^l 
between MHD and HD cases. This is a measure of the reduced horizontal mixing of internal 
energy between up- and downflows. We find a significant reduction of the horizontal energy 
mixing with the presence of strong magnetic field. This results in an increased and more 
efficient energy transport, i.e. the enthalpy flux is not suppressed significantly, despite the 
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reduction of overturning mass. 


4. Summary 


We presented a series of high resolution MHD simulations in the solar convection zone 
and studied small-scale dynamo action through the convection zone. For grid spacings 
smaller than 350 km, the RMS magnetic field around the base of the convection zone reaches 
95% of the equipartition field strength. Even around the top boundary 55% of equi partit ion 


magne tic field is maintained. These values are consistent with those required in 


Re mpel 


(120141) to explain the observed field strength of the quiet Sun. We find a significant Lorentz 


force feedback, the convective RMS velocity is reduced to 50% and 85% compared to the 
hydrodynamic reference run at the top and bottom boundaries, respectively. 


The efficiency of the small-scale dynamo achieved in this study is larger than that of 
the large-scale dynamo found in many recent global simulations (i.e. we reach substantially 
stronger B rms values in our setup with a small-scale dynamo alone). 

The significant suppression of the convective velocity does not lead to a significant 
decrease of the enthalpy flux. This study shows that the suppression of the velocity reduces 
the horizontal mixing of the entropy between upflow and downflow regions and causes larger 
entropy perturbations. This improves the efficiency of the convective energy transport. The 
reduction of the kinetic energy does not converge especially in the lower convection zone 
where the convective scale is larger, which implies that Lorentz force feedback on convective 
flows is potentially larger than presented here. Even in that case the horizontal mixing 
of the entropy might be further reduced leading to even stronger entropy perturbations 
between up- and downflows and more efficient convective energy transport. Overall we 
conclude that the Lorentz force feedback from small-scale magnetic field is likely significant 
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and has to be taken into account for realistic simulations of the solar convection zone, which 
is currently not possible in global simulations. 


Recently 


Hanasoge et a l. 


(12012 ) presented rather stringent helioseismic constraints on 


convective velocities in the solar interior, which suggest that current convection simulations 


overestimate convective velocities by up to several orders of magnitude. In a co mplete! 


independent study of solar supergranulation through numerical simulations 


Lord et al 


y 


(2 0141) found that high resolution simulations of the solar photosphere are consistent with 
observational constraints in the uppermost 10 Mm of the convection zone, but overestimate 
convective velocities by a factor of a few in more than 10 Mm depth. The work presented 
here indicates that ubiquitous small-scale magnetic field throughout the solar convection 
zone might contribute to some degree to resolving the issue, however, the reduction of 
convective velocities we find is still insufficient. Not completely unrelated to this is likely 
also the difficulty global convection simulations have in reproducing a solar-like diffe rentia l 


rotation in set ups that transport a full solar luminosity. The work by 


Gastine et al. 


Kapvla et al. 


(2014j); 


( 20141 ) suggests that a fast rotating equator requires a lower Rossby number, 


i.fu lower convection velocities throughout the convection zone. Recently 


Fan fe Fang 


(12014 ) showed that also Maxwell-stresses from a dynamo generated magnetic field can 


significantly alter the sign and shape of differential rotation. To which degree this is due 
to large-scale or small-scale magnetic field components remains at this point an open issue 
and requires numerical simulations of efficient small-scale dynamos in a global setup with 
rotation, which is work in progress. 
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Fig. 1.— Contours of (a) vertical velocity v x in the unit of m s 1 and (b) vertical magnetic 
field B x in the unit of G are shown at r = 0.95 R@ for case M256D. 


Table 1: Number of grid points and diffusivities in each case. 


Case 

N x x N y x N z 

k, v, r) [cm 2 s *] 

H256D, M256D 

72 x 256 x 256 

1 x 10 12 

H512, M512 

144 x 512 x 512 

- 

H1024, M1024 

288 x 1024 x 1024 

- 

H2048, M2048 

576 x 2048 x 2048 

- 
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Fig. 2.— Temporal evolution of magnetic energy averaged over the whole calculation domain 


in case M256D. 
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(a) 



X/Rsun 


(b) 



X/Rsun 


Fig. 3.— (a) Black lines shows the equipartition magnetic held to the kinetic energy ( B cq = 
V^rrms). Solid and dotted lines shows results in MHD (M256D) and hydrodynamic 
(H256D) cases. Red line shows the RMS magnetic held, (b) The ratio of R rms to B cq in case 
M256D. 
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Fig. 4.— The distribution of vertical velocity v x in case H2048 at r = 0.95_R o (panels a and 
b) and r = 0.8i? e (panels c and d) in the unit of m s -1 . Panels b and d show the zoom in 
image of panels a and c, respectively. The white rectangles indicate the region for panels b 
and d. 
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Fig. 5.— Spectra of kinetic energy at (a) r = 0.95 R 0 , (b) r = 0.8 R e , and (c) r = 
0.72 Rq. Green, black, blue, red lines show the results of H256D, H512, H1024, and H2048, 
respectively. The dotted lines show a Kolmogorov slope of -5/3. 
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Fig. 6.— Temporal evolution of magnetic energy. The black, blue, and red lines show the 
results in cases M512, M1024, and M2048. 
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Fig. 7.— Contour of vertical velocity v x (panels a and b) in the unit of m s _1 and vertical 
magnetic field (panels c and d) in the unit of G are shown at r = O.95_R 0 for case M2048. 
Panels b and d are the zoom-in figure for the panels a and c, respectively. 
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Fig. 8.— The same format as Fig. [7] at r = 0.8 Rq. 
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Fig. 9.— Joint PDFs of vertical velocity (y x ) in the unit of m s 1 and vertical magnetic 
field (. B x ) in the unit of G in case M2048 at (a) r = 0.95-Rq and (b) r = O.8i? 0 . 
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Fig. 10.— Spectra in cases H2048 and M2048. Panels a, b, and c denote the result at 
r = 0.95 Rq, 0.8 R, and O.72i? 0 . The black and blue lines show kinetic energy, and magnetic 
energy in case M2048. The red line shows the kinetic energy in case H2048. Panel d shows 
the ratio of E mag to E kin at r = 0.95 R e (black), 0.8_R e (blue), and 0.72 R & (red). 
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(a) 



(b) 



Fig. 11.— The relation between equipartition magnetic field B eq and RMS magnetic field 
B rms . (a) The solid and dash-dotted lines show B cq and R rms , respectively for the cases M512 
(black), M1024 (blue), and M2048 (red). The dotted black, blue, and red lines show B cq in 
cases H512, H1024, and H2048, respectively, (b) The ratio of R rms to B cq in cases M512, 
M1024, and M2048. 
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Fig. 12.— The work densities of pressure/buoyancy {W&. solid line) and the Lorentz force 
(W\. dotted line) (see eqs. (IT2]) and (fill) ). The black, blue and red lines show the results of 
M512, M1024, and M2048, respectively. 
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(a) 


(b) 





Fig. 13.— The solid (dotted) black, blue and red lines show the results in cases M512(H512), 
M1024(H1024), and M2048(H2048), respectively, (a) The distribution of RMS velocities, (b) 
The ratio of total RMS velocity urms in the MHD cases to hydrodynamic cases, (c) The ratio 
of vertical RMS velocity w^rms in the MHD cases to hydrodynamic cases, and (d) The ratio 
of horizontal RMS velocity UhRMS { v h = + v%) in MHD cases to hydrodynamic cases. 
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Fig. 14.— Horizontally integrated unsigned mass flux po\v x \. The solid (dotted) black, blue 
and red lines show the results in cases M512(H512), M1024(H1024), and M2048(H2048), 
respectively. 
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Fig. 15.— PDFs for (a) the vertical velocity v x , (b) a horizontal velocity v y , (c) the vertical 
magnetic held B x and (d) a horizontal magnetic held B y . The black, blue and red lines show 
the result at r = O.95i? 0 , O.8i? 0 , and O.72i? 0 , respectively. The dotted and solid show the 
result in cases H2048, and M2048. 








40 


(a) 



x/R® 


(b) 



x/R® 


Fig. 16.— (a) Kurtosis and (b) skewness of the velocities. Black, blue and red lines show 
v x , v y: and v z , respectively. Dotted and solid lines show the results of H2048 and M2048, 
respectively. 
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Fig. 17.— Energy fluxes in cases H2048 and M2048 (dotted and solid lines, respectively). 
The black, blue, red, light green, and purple lines show the total, kinetic, enthalpy, radiative, 
and Poynting fluxes, respectively. Black solid line shows the F — 0 axis for a reference. The 
solar flux in this study is F Q = 1.21 x 10 11 erg s 1 cm 2 . 
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Fig. 18.— PDFs of the entropy. The black, blue and red lines show the result at r = 0.95 -Rq, 
0.8-Rq, and O.72i? 0 . The dotted and solid lines show the results in cases H2048 and M2048, 
respectively. 
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(a) 



(b) 



Fig. 19.— Joint PDFs with the vertical velocity v x and the entropy si at r = 0.80-R©. The 
panels a and b shows the result in cases H2048 and M2048, respectively. 
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Fig. 20.— Contour of entropy at r = 0.8-R© in cases (a) H2048 and (b) M2048 in the unit 
of erg g -1 K _1 
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Fig. 21.— Spectra of entropy at (a) r = 0.95 R Q , (b) r = 0.8 R Q , and (c) r = 0.72i?©. The 
red and black lines show the result in cases H2048 and M2048, respectively. 




Fig. 22.— Spectra of enthalpy flux at (a) r = 0.95 R Q , (b) r = 0.8 R e , and (c) r = 0.72 R & . 
The red and black lines show the result in cases H2048 and M2048, respectively. 
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(a) 


(b) 




Fig. 23.— (a) The RMS values of entropy perturbation. Black, red, and blue lines show 
the mean value, value in upflow and downflow, respectively. Dotted and solid lines show 
the results in cases H2048 and M2048, respectively, (b) The ratio of the RMS entropy 
perturbation between the cases H2048 and M2048. 
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Fig. 24.— (a) Kurtosis and (b) Skewness of the entropy. Dotted and solid lines show the 
result in cases H2048 and M2048, respectively. 
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Fig. 25.— Ratio of horizontally averaged |siv y | between MHD and HD cases. The black, 
blue and red lines show the result in cases H(M)512, H(M)1024, and H(M)2048, respectively. 








